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Abstract 

We have performed a many-body calculation of the longitudinal nuclear (e, e') re- 
action employing a Second RPA (SRPA) formalism which contains the A(1232). More 
explicitly, our scheme contains RPA correlations as well as Hartree-Fock and second 
order self-energies, where an accurate evaluation of exchange terms is achieved. Using 
this formalism we have evaluated the longitudinal response function for 40 Ca. We give 
final results at momentum transfers ranging from 300 up to 500 MeV/c, obtaining a 
good agreement with data. 
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1 INTRODUCTION 



The quasi-elastic peak in inclusive electron scattering for medium and heavy nuclei has been 
paid much attention in recent years. The good description of the cross section in the former 
work of Monitz et al. [1] was suddenly broken when the longitudinal and transverse responses 
were experimentally separated by Meziani et al. [2]. Even though the total cross section is 
well reproduced by a Fermi gas, this simple model overestimate data in the longitudinal 
channel, while it does the opposite in the transverse one. More recently Williamson et al. [3] 
have presented new measurements with an increase (decrease) of the longitudinal (transverse) 
nuclear response. Still, the discrepancy remains. 

Many theoretical efforts have been developed to overcome this problem. Let us resume 
the main ones. One group ascribes the problem to modifications of the nucleon properties 
within the nuclear medium by the employment of small effective masses, swelling of nucleons, 
abnormal nucleon form factors, etc. (see Ref. [4]). Even though this approach is successful 
for the longitudinal channel, it fails in the transverse one. The other group, to which belongs 
the present work, is based on the many-body theory. 

Many-body calculations were performed in finite nucleus [5]- [14] or in nuclear matter 
[15]- [30]. Fortunately, surface effects for heavy and medium nuclei in the energy-momentum 
region of interest are not very important and the nuclear matter formalism is a good ap- 
proximation for the nuclear system, once a variable Fermi momentum or the local density 
approximation is used. On the other side, the many-body problem is very complex: ini- 
tial and final state interactions, meson exchange currents, etc., should be considered, which 
implies a huge numerical effort, even in nuclear matter. 

Perhaps the simplest way of introducing nuclear correlations is the RPA (see Refs. [6], 
[11], [20] and [21]). In the same line, a step further in difficulty is the SRPA [7], [24], which 
opens decay channels beyond one particle-one hole (lplh) states and the Extended RPA 
(ERPA) ([8], [19] and [22]), adding more complex ground state correlations than in the 
RPA. Others approaches are the Correlated Basis Function (CBF) [9], the Green function 
method [10] and the one Boson Loop Expansion (BLE) [29], [30]. Looking at these formalism 
in terms of Feyman diagrams, all these approaches should converge to the same result, once 
the proper number of diagrams are considered. In spite of the sophistication of all these 
models, a clear understanding of the nuclear response at any momentum transfer is still 
elusive. 
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In some recent works the role of the A(1232) on the longitudinal channel was explored in 
some recent works. Gil et al. [28] study both channels within an elaborate formalism which 
includes two-body corrections to the excitation operator, initial and final state interaction 
and the A. The A was considered by modifications in the external operator, self-energy 
insertions and in the nuclear interaction by means of the ring series with the A (polarization 
effect). The A is also incorporated in the longitudinal response in the works of Cenni et 
al. [29] and Amore et al. [30] where the Boson Loop Expansion is employed. In Amaro et 
al. [18] (and also in Ref. [30]), it is addressed that a A — hole pair can be directly excited in 
the longitudinal channel. In Ref. [21] we have also examined both longitudinal and transverse 
responses with the explicit inclusion of the A. In that work, however, we concentrate on the 
transverse channel, with special attention on the RPA-exchange terms, concluding that an 
accurate evaluation of them is important. In that early work, the A has entered into the 
longitudinal response only through the second order self-energy, as there were no RPA terms 
in the longitudinal channel due to the particular election of the residual interaction. 

In the present work we have developed a SRPA formalism which employs a more elaborate 
residual interaction than in our previous work [21]. The SRPA contains the RPA correlations 
together with the Hartree-Fock and second order self-energies. Exchange terms were carefully 
taken into account and the A was incorporated both in the many-body calculation and by 
it direct excitation. 

The work is organized as follows. In Section II, we present the formalism. In Section III 
we analyze the effect of each ingredient in our model and analyze results for m Ca. Finally, 
in Section IV some conclusions are given. 
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2 FORMALISM 



Even though the main ingredients of our formalism were already presented in Refs. [21] and 
[24], here we have preferred to show a complete version of it. This is because there are some 
elements which are particular to the nuclear longitudinal response with the A which should 
be discussed. 

The nuclear response function to an electromagnetic probe in the longitudinal channel is 
defined as, 

R L (q, u) = --Im < \0[G(lu)O l \ >, (1) 

7T 

where q represents the magnitude of the three momentum transfer by the electromagnetic 
probe Ol, uj is the excitation energy and | > is the uncorrelated nuclear ground state. 
Ground state correlations beyond RPA are not analyzed in this work. The polarization 
propagator is given by, 

GH = u-H + ir, " u + H-ir, ' (2) 
where H is the nuclear Hamiltonian. As usual, H is separated into a one-body part, H , 
and a residual interaction V. 

We present now two projection operators P and Q. The action of P is to project into 
the ground state, the one particle-one hole (ph) and one delta-one hole (Ah) configura- 
tions. While Q projects into the residual n p particle-n^ hole-nA delta configurations. More 
explicitly, 

P=\X\+P N + P A , (3) 



with 



P N ="£\lplh><lplh\, (4) 



lplh 



P A = \ lAlh >< lAlh \ ( 5 ) 



and 



Q = \n p p (n h - n p )A n h h >< n p p (n h - n p )A n h h\, (6) 

n h > 2 
< n p < rih 

where we have introduced P/v and Pa for convenience. It is easy to verify that P + Q = 1, 
P 2 = p, Q 2 = Q, and PQ = QP = and also, P i P j = 5^ and PQ = QP i = (i = N, A). 
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We insert now the identity into eq. (1). Note that the external operator can connect the 
uncorrelated ground state only to the P space if we use a one body operator for Ol, we 
have, 

R L (q,uj) = --Im< \0[ P G pp (lo) P O l \ >, (7) 



7T 

where Gpp = PGP. It is easy to see that, 

1 1 



Gp P (u) 



u - H PP - £ p Q p + ir] co + H PP + £ p Q p -ir]' w 
where, 

S PQP = Vpn Vqp - Vpn ^ Vqp, (9) 

with obvious definitions for Hpp, etc. 

Until now, the formulation is quite general. The next step is to adopt a model for the 
external operator, the nuclear Hamiltonian and the nuclear states. The nuclear interaction 
is discussed in the next section, while the nuclear structure problem is the main issue of this 
section. Before this, we show a model for Ol, which is divided into two contributions, 

o L = o ph + o Ah , (10) 

where the action of O ph (0 Ah ) is to create a ph (Ah) pair. The matrix elements for O ph and 
Oj\h were taken from Refs. [17] and [18], respectively, 



< lplh\O ph \ >= G p E (q,u)^=xi p Xt p ^^Xs h Xt h 

and 



(11) 



< lAlh\0 Ah \ >= G A (q,cu)^^xt A xLS* ■ (h x q) T 3 XshXth (12) 



m 2 

where the s's and the t's stand for the spin and isospin quantum numbers, h is the momentum 
carried by the hole, k = q/(2m), r = (q 2 - uj 2 ) / (Am 2 ) , r' = ^ (r + ( m 2 A mA m ) 2 ), and m (m A )is 
the nucleonic (A) mass. In eq. (12) the Pauli matrices cr and r 3 were replaced by the 
corresponding transitions matrices S and T 3 [31]. The electromagnetic form factors are, 

and 

G AM= 2M^ M (l-^^T^ (14) 
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In the longitudinal channel there is no interference term between O p h and A h- This is at 
variance with the case of the transversal one and it is a consequence of the presence of the 
hole momentum in A h- This property allows to rewrite Eq. (7) into two terms, 

Rl = R P h + Rah, (15) 

where 

R ph = --Im< \0\ h P N Gp nPn {uo) P N O ph \ > (16) 

71 

and 

Rah = -- Im< \Oi h P A G PaPa (lu) P a Ah \ > . (17) 

71 

In the next two subsections we analyze separately R p h and R A h- 
2.1 R p h contribution 

The starting point of this subsection is Eq. (16) where our concern is to understand this 
expression in terms of it physical ingredients. To start with, we replace Pjv by hs explicit 
expression, 

= ~\lm £ < |0> P lft X iplftl _ x _ llp'iy > 

" lplh,lp>lh> ^ ^PnPn ^ T ill 

x < lp'lh'\O ph \ > (18) 

where for simplicity we have shown only the forward going contribution. We consider now 
the energy denominator of Eq. (18). For H PnPn we have, 

< lplh\H PNPN \lp'lh' >= 6 ph , p , h ,(Ei% + Eg* ) + V lplh , lp , lh ,, (19) 

where E[° p \ h is the kinetic energy of the particle-hole pair, is the Hartree-Fock self- 

energy contribution and Vi p ih,i p >ih' is the residual interaction between two particle- holes 
pairs. More explicitly, we can define the energy of a lplh pair up to first order in the 
residual interaction as, 

<k-< ) i, + Sfi = e«(p) - e «(*0, (20) 
where in addition, we have introduced the energy of a single nucleon with momentum p as, 

= +E* F (P). (21) 



The analytical expression for H HF {p) depends on the interaction employed and it is discussed 
in the next section. 

The action of the T i PnQPn self-energy is to connect the Pjy-space with the Q one. This 
means that £ PjvQPjv opens the decay channel of a lplh pair into a 2p2h configuration, a 
lAlp2h one and also more complex configurations like 2A2h. 

The matrix elements of £ PjvQPjv are written as, 

< l P lh\V™ p »\lpflh' >= 5 ph , plhl ^ PN (h,q,u;) (22) 

where we have kept only the diagonal terms of the self energy. This approximation is based on 
numerical reasons already discussed in Ref. [21]. Also, in this work we have studied in detail 
Sf^vT ^ splitting it into four components, depending on the character of the intermediate 
configuration Q, as, 

T^{h,q, U ) =j:^S iPN (h,q,u). (23) 
i=i 

where Qi represents a lplh bubble attached to a lplh, Q 2 is a lAl/i bubble attached to a 
lplh, Qs is a lplh bubble attached to a lAl/i and finally, Q4 is a lAlh bubble attached to 
a lAl/i. Diagrams of each of the Qi can be found in Fig. 6 of Ref. [21]. Their analytical 
expressions are also found in the same reference. 

Also through the numerical analysis it turns out that the dependence of the self-energy 
over the hole momentum is not very strong. This allows us to make an average over it, 

^ Q h PN (h,q,u;) = T$g iP »(q,u>) (24) 

where 

X%& iP »{q,u) = fd 3 h 9(k F - \h\) ^f P -{h,q,u) (25) 
3 

and k,F is the Fermi momentum. Now Eq. (18) reads, 



Rn 



-- Im ]T < \°lh\ l P lh >< l P lh \ 



vh = nn 2^ ^ ^ ipi"-\ " =p) 

71 lplh,lp'lh' W — Oph,p'h' ^\p\h ~ V\p\h,\p'\h' + 

x\lp'lti >< lp'lh'\O ph \ > (26) 



where 

Spl/i = ^lplh 1" ^lplh 1" ^lplh 



^lplh - ^lvlh + h lplh + ^lplh {Q,V). {2() 
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The non-diagonal term of the interaction Vi p ih,i p 'ih' , is what makes Eq. (26) hard to evaluate. 
Due to this reason, it is convenient to show first the result without Vipifc,ip'ih' ■ We define, 

R pI = -^ Im J2 I < Wh\O ph \ > | 2 -i (28) 

where the lplh energies E^ lh with % — 0, 1 and 2, were already defined in Eqs. (19, 21) and 
(27), respectively. For i = 0we obtain the free response, 



<W) = E I < ^h\O vh \ > I 2 6(u, - (^1! _ (29) 

lplh P " 



where h + q is the particle momentum and we have allowed the particles and the holes to 
have different masses. Performing the summation over spin and isospin and making the 
conversion of sums over momenta into integrals we have, 



Rph(Q,u) = {G E {q,u)) 2 ^- £{q,uj,m p ,m h ) (30) 



where A is the mass number and C(q, ou, m p , m^) is a modified Lindhard function defined in 
Appendix A. The Goldstone diagram of R p ^ is drawn as the first one in the r.h.s. of Fig. 1. 
The response function for i — 1 is, 

= E I < Wh\O ph \ > | 2 S(u - (e^(\h + q\)- e W(h))) (31) 

lplh 

The expression within the 5-function is a transcendental equation due to the presence of 
Tip/f and is solved numerically. Alternatively, in Ref. [27] it was discussed that the 
response is well reproduced by R p ^ but using two effective masses (one for particles and one 
for holes), and an energy shift. Adopting this approximation, we can write, 



Rph\<l,u) = (G E (q,u)f ^- C(q,uj,m p ,m h ) (32) 

where u = uj + Acj and the values for m p , rrih and Au depend on the interaction and on the 
momentum transfer. The first three diagrams contributing to R^ are the first three ones of 
Fig. 1. The fact that the Hartree-Fock self-energy is in the energy denominator implies that 
this contribution is sum up to infinite order. 

(2) (1) 

Finally we consider R K £. At variance with the case of R p ^ where the the Hartree-Fock 
self-energy is a pure real function of the particle (or hole) momentum, the T^^ QPn self-energy 
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has a real and an imaginary part. It is easy to see that R p 2 ^ can be written as, 

E m (a rJFR^(aF) /mgSg^oj 

(33) 

The lower order contributions which comes entirely from E^^* Pjv , are shown in the second 
line of Fig. 1 (in upward order of i). 

We turn now to the discussion of R p h with the inclusion of the V\ p \h,\p'\h' term. We write 
down from Eq. (26) the matrix element of the polarization propagator, 

G(u) lplh}lp > lh , = — — — (34) 

u — Oph,p'h' ^ipih ~~ Vi p ih,ip>ih' + iV 

where as mentioned above, we show for simplicity only the forward going contribution. To 
treat this, the standard Dyson equation is employed 

G{u})i p \h,ipfih' — G( 2 \uj)ipih + G^^ipihVipih^piih'G^ipi^ipiih' , (35) 
where we have used G^ 2 \uo)i p i h defined as, 

G (2 \u) lplh = -1 (36) 

^ - E lplh + w) 

instead of the free polarization propagator G^(uo) ph . By successive iterations Eq. (35) leads 
to a series which can be summed up only in some few particular cases: when exchange 
terms are neglected or when V\ p \h,\p'\h' is a contact or a separable interaction. In the case 
that exchange terms are neglected the series is known as ring series and the sum of it as ring 
function. The puzzle is then, how to deal with exchange terms. For a general force, exchange 
terms should be evaluated numerically order by order and in practice, this can be done up to 
second order in the residual interaction. Third and higher order contributions, imply a quite 
involve numerical task due to the number of diagrams and the dimensionality of the integrals 
required to evaluate each diagram. However, exchange terms of a contact interaction can be 
easily included by a re-definition of the constants entering into the interaction. Using these 
facts, in Ref. [20] we tackle this problem as follows: we re-write the residual interaction by 
summing and subtracting a contact force, V\ p \h,\p>\h' ■> 

V\ p ih,ip'ih' = Vipih,ip>ih> + ^ipih,ip'ih' (37) 
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where Vf plh 

lp'W = Vipih,ip'ih' — Vv p ih,ip'ih' contains the finite range part of the interaction. 
Varying V^ lhlp , lh ,, the interaction V^ lh lp , lh , is adjusted in order to ensure a fast convergence 
of the exchange terms. More specifically, the different constant terms of Vf plhlp , lh , (that is, 
the ones of V^ lhlp i lh i l ) , are varied until the second order exchange contributions are almost 
negligible compared with the first order ones. 

The polarization propagator of Eq. (35) can now be drawn as 

G(u)i plhtlp >i h > = G(u)i plh + G(uj)( plhlp , lh , + G((j)i p F lhlp , lhl , (38) 



where, 



r C _ r (2) ,^(2) V C r (2) , r {2) V C r (2) V C r (2) 

^lplh — ^lplh ^ lplh v Iplh^lplh^ Iplh (j lplh v lplh,lplh (j lplh v lplh,lplh^ r lplh 



(39) 



r F _ r (2) T/ F r {2) r (2) T/ F r {2) tt-F r (2) / m \ 

^lplh,lp'lh' — { ~ sr lplh v lplh,lp'lh' { ~ :, lp'lh' "r ^ Jr lplh v lplh,lp"lh" (j lp"lh" v lp"lh" ,lp'lh ,(j lp'lh' ■ l^ u J 



r CF _ r (2) rrC r (2) i/F r (2) , 

~~ [ ~ r lplh v lplh,lp"lh" { ~ r lp"lh" v lp"lh",lp'lh' { ~ r lp'lh' "T 

r (2) T/ F r (2) V C r {2) , 

V lp"lh",lp'lh l(j lp'lh' "i" 

, ^(2) T/ C ^(2) T/ C ^(2) T/ F ^(2) 

^lpl/i V lplh,lp'"lh" l(j lp'"lh'" V lp"'lh"',lp"lh" (j lp"lh" V lp"lh",lp'lh' i ~ r lp'lh' 



(41) 



Inserting Eq. (38) into Eq. (26) one can define three different contributions to the response 
function, R$ h , R F h and i?gf, associated to Gf plh , Gf plhAp , lh , and Gg£ Mjm ,, respectively. 
In the case of R ph , since V-^pih^ipiih' is a pure contact interaction, the direct and exchange 
terms are equivalent and they can both be summed up to infinite order by evaluating the 
ring series of Eq. (39). In the case of R ph) Eq. (40) is only considered up to second order 
in V-^ lhlp , lh ,, while for R p F , Eq. (41) is also evaluated up to infinite order in the V^ lhlp/lh , 
interaction keeping terms up to second order in Vf plh lp , lh ,. Our final result is then, 

Rph = Rp h + Rp h + Rph (42) 



Eventually, Vi p ih,ip'ih' can contain constant terms, but these terms can not be changed. Only the ones 
of Vi plh lp , lh , can be varied, as this interaction is an artificial device to guarantee the fast convergence of 
exchange terms evaluated with Vf plh lp , lh , 
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Explicit expressions for the different exchange terms needed to build up the R p h contribution 
are very similar to the corresponding ones in the transverse channels reported in Ref. [20]. 
There are minor differences in the multiplying constants, the spin and isospin sums and the 
energy E[°J lh should be replaced by Ef p \ h . For this reason, we do not reproduce them. We 
neither show the ring series nor the direct terms as they can be found in many references 
(see for example Ref. [32]). If in Eq. (42) we replace G^(uj) lplh by the free propagator 
ipih, the RPA approximation is obtained. The first order direct contribution to the 
RPA is drawn as the fourth diagram of the first line in Fig. 1, while the next diagram is 
the corresponding exchange term. Our approximation, where self-energies are included in 
G^(u) p h, is usually known as Second RPA (SRPA). 

2.2 R^h contribution 

At variance with R ph) the particular character of A h does not allow a lAl/i bubble to 
propagate as a RPA-series. In this case, only the self-energy modifies the free response. The 
expression for R Ah is obtained following closely the steps done in the last sub-section to 
arrive to R^- Obviously, one should replace Pv by Pa an d do some others minor changes. 
Before presenting the expression for R Ah , we would like to discuss briefly the lower order 
contributions to the Raii response shown in Fig. 2: the first diagram in the r.h.s. of the 
first line is the free response. The second and third diagrams represent the Hartree-Fock 
self-energy over a hole line. There is no Hartree term for the A due to the absence of a 
scalar-isoscalar term in the ViAife,iA'ih' interaction. The fourth diagram in the first line 
of this figure is the Fock contribution to the A. Finally, in the second line, we show the 
diagrams with self-energy insertions ^i A ih A ■ 

Now we briefly outline the main expressions required to obtain R^ Ah - In Eq. (17) we 
substitute Pa by it explicit expression, 

R Ah = --Im <|OLJ1A1/i><1A1/i| 77 VpaQPa-l. • \ lA ' lh '> 

x < lA'lh'\0 Ah \ > (43) 

The energy denominator of Eq. (43) is approximated as, 

< lAlh\H PAPA + Z p ^ p *\lA'lh' >^ 5 Ah , Alhl (E[°l lh + Ef/ lfc + Sf^ p -). (44) 
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It is convenient to define the energies, 

E?L = E^ + Z^^u,) (45) 

where expression for is given in Appendix B. We have adopted the same approxima- 

tions for SfA^T A as for S^^ Pjv . Explicit expressions for both self-energies can be found in 
Ref. [21]. 

From Eq. (43) we define now, 



R% = --ImJ2 | < lAlh\0 Ah \ > | 2 -i (46) 



After some algebra it is easy to find, 
,(o) , x 8n 2 A kp 



RZ(q, uj) = (G A (q, u)f (1 + r') C A (q, u, m A , m h ) (47) 



B£L(q, u) = 8n2 ^ 4 kF (G A (q, uo)f (1 + t') C A (q, Q, m A , m h ) (48) 



ffP) ffl ^_ 1 r dFR W( aF ) ImX™ p *(q,u>) 

(49) 

where in analogy with the p/i-case, we have employed a modified function £ A (q, cu, m A , nth), 
which is defined in Appendix A. 
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3 RESULTS 



In this section we show explicit results for the nuclear quasi-elastic longitudinal response. The 
properties of medium mass nuclei in the energy-momentum region of interest are reasonably 
described by non-relativistic nuclear matter once a proper Fermi momentum is used [16]. In 
particular, we consider the response of 40 Ca using a Fermi momentum hp = 235 MeV/c. 

For the residual interaction V\ p ih,ip>ih' ■> we have made use of the parameterization of the 
Bonn potential [33], as presented in Ref. [27]: it is expressed as the exchange of n, p, a and 
uj mesons, while the r\ and 5-mesons are neglected. Also from this reference, we have taken 
the expressions for the Hartree-Fock self-energy When at least one A is involved, 

we have adopted the model interaction already used in Ref. [21], which for convenience is 
reproduced in Appendix B. Employing these interactions, the first step is to solve Eq. (26) 
(and the corresponding one in the A sector, R^ h ), in order to fix the effective masses and 
energy shifts for particles, holes and deltas. The results for several momentum transfers is 
shown in Table I. 

Our main concern is to understand the longitudinal response and in particular the effect of 
the A over it. As already stated in the last section, the A affects the longitudinal response via 
two mechanisms: by the direct excitation of a A/i-pair and through the self-energies Sf^^ Pjv 
and Sf^^ PA with i — 2, 3 and 4. The first effect is considered in the R^h response function. 
In Fig. 3 we present results for R^ h , R^ h and R^ h at momentum transfer q = 400 MeV/c. 
The comparison of R^ h with R^ h shows that the effect of the Hartree-Fock self-energies is 
not relevant. At variance, the influence of the S p ^^ Pa is very significant: it produces an 
important redistribution of the intensity to lower energies. The S p ^^ Pa self-energy has a 
real and an imaginary part and both terms are connected through a dispersion relation. The 
effect of the imaginary part is to spread the response, while the real part shifts the energy 
position of the peak to lower values. These effects are very strong because S p ^^ Pa increases 
with energy and the Ah peak is high in energy. It will be shown soon that when we compare 
the magnitude of RAh with R p h and with data, the R^h response is as a whole negligible. 
However, R& h had been discussed for completeness. 

In Fig. 4 we analyze the effect of the different contributions needed to built up the R ph 
response. In each panel, we have drawn the free response, R p °h, the total response R p h, and 
the experimental points. This was done as a guidance to understand the behaviour of each 
term separately. Let us resume the meaning of dashed lines for each panel: 
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Panel (a): dashed lines represents the response. Which means to study the effect of 
the Hartree-Fock self-energy. 

Panel (b): we show the RPA response; that is, in Eq. (42) we have employed G^\u) p h 
instead of G^ 2 \u) ph . 

Panel (c): here we isolate the effect of the S0^^ lPjv self-energy. In this panel, we have 
evaluated Eq. (33) using a re-defined energy, E^ lh = Ef p \ h + Eg^ lPjv (q, uj). 

Panel (d): we evaluate Eq. (42) excluding the A terms of the self-energy E^^ Pjv (which 
means to use Sf^ lPjv alone). 

The Hartree-Fock self-energy spread the free response and moves the quasi-elastic 

peak to higher energies, as shown in Panel (a). This result was already reported in others 
works (see for example Ref. [27]). The effect of the RPA (see Panel (b)), depends on the 
character of the interaction. In the longitudinal channel, the direct contribution to the RPA 
(ring series), selects the isospin scalar end isospin vector central terms of the interaction 
(usually named as F and F'). For our election of the interaction, these terms are repulsive, 
which is reflected in the response. Turning now to Panel (c), we see that the action of the 
^ipife lPjv self-energy is similar to the one in the A-sector. However, the magnitude of the 
spread and energy shift are smaller. Perhaps, the most interesting result of this figure is the 
one in Panel (d). The comparison of our final result with the dashed line, indicates that the 
A is very important in the longitudinal channel, not because of it direct excitation, but due 
to its indirect effect through the self-energy. 

In Fig. 5 we compare our final result with data for several momentum transfer q, rang- 
ing from 300 up to 500 MeV/c. We have obtained a good agreement at low values of q. 
Although the agreement for q = 450 and 500 MeV/c is not so good, the accordance with 
the total intensity is satisfactory. At this point it is important to mention that the values 
of the parameters entering into the residual interaction when the A is involved are quite 
uncertain. This problem was already discussed in Ref. [21] where we showed that self-energy 
contributions which contains the A, change considerably due to this ambiguity. Eventually, 
and in view of the result of Panel (d) in Fig. 4, a more accurate interaction could improve the 
agreement for high values of q. In addition, the p/i-interaction itself is not unambiguous. Let 
us quote Ref. [27] where it has been pointed out that the bare nucleon-nucleon interaction 
employed is too repulsive to reproduce the correct binding energy and this shortcoming might 
also affect the RPA spectrum (overestimating the hardening of the response). In summary, 
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although our results are encouraging and in fact we have obtained perhaps one of the best 
agreements for the longitudinal response in a wide region of q; these results are conditioned 
by the election of the residual interaction which is yet unknown. 

Concerning a comparison with other approaches, it is interesting to comment on two 
works already mentioned in the Introduction, where the longitudinal response with the in- 
clusion of the A is considered. The first one is due to Gil et al. (Ref. [28]), where both the 
longitudinal and transverse responses are evaluated. This work is very complex in the sense 
that most of the many-body mechanisms where incorporated: mesons exchange currens, po- 
larization, initial and final state interactions and the A. The comparison with data is done 
only at two momentum transfers: q = 300 MeV/c for 12 C and q = 410 MeV/c for A0 Ca. 
The last results are shown for energies ranging between 40 and 220 MeV. Even if there is 
a small tendency to overestimate data, these results certainly represents an improvement 
over previous ones. Unfortunately, final values for 40 C a are shown for only one momentum 
transfer and in a narrow energy region. In addition, the role of the A in the longitudinal 
channel is not discussed in particular. The second work was done by Amore et al. (Ref. [30]), 
where they have studied the effect of the A over the longitudinal response using the BLE. 
As in our case, they have considered both the direct excitation of a A/i-pair and the effect 
of the A over the response through many-body contributions. They have also incorporated 
relativistic effects and done the comparison with 12 C data for q = 400, 500 and 600 MeV/c. 
Their final results underestimate the experimental points. Using a different formalism and 
interaction, this work agrees with our results in two main issues: the direct excitation of a 
A/i-pair is not significant and the many-body effect of the A is important. 

The main difference between our scheme and those of Refs. [28] and [30] is the way in 
which we deal with the Pauli exchange terms. In several works (see for instance Refs. [21], 
[23]), we have determined that exchange terms are important and need to be accurately 
evaluated. If we now incorporate the transverse channel into the discussion, there are many 
contributions like meson exchange currents and ground state correlations beyond the RPA, 
which also need to be considered. However, in the longitudinal channel these terms are less 
relevant than in the transverse one. Moreover, the interaction utilized in Refs. [28] and [30] 
is different than ours (in fact, each work employs a different one). This point is particularly 
important not only because of the sensibility of the results over the interaction, but also 
because antisymmetrization establishes constrains over the interaction (see Ref. [34]). Let 
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us be more specific about these constrains. When an antisymmetric matrix element is 
evaluated, the antisymmetrization operator can act either over the initial and final state or 
over the interaction itself. It is possible to build up an 'antisymmetric' interaction where 
the prescription is to evaluate direct matrix elements with this force. This is equivalent 
to employ the non-antisymetrized interaction with direct and exchange matrix elements. 
When the antisymmetric interaction is employed, the different channels of the interaction 
(for instance, F, F', G and G' in a Landau-Migdal type interaction), are not independent: 
a change in one affects all the others. As the RPA longitudinal and transverse responses 
are dominated by different channels of the interaction, the just mentioned constrains should 
be carefully considered when dealing with both channels. Alternatively, the employment 
of any non-antisymetrized interaction within a scheme which accurately evaluates exchange 
terms is equivalent, as just mentioned. In a forthcoming work we will discuss this issue in 
detail [26]. 
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4 CONCLUSIONS 



In the present work we have analyzed the longitudinal response of the inclusive electron scat- 
tering process. We have presented a formalism developed for non-relativistic nuclear matter, 
which includes correlations of the SRPA type. The A-degree of freedom was incorporated 
via two mechanisms: the direct excitation of a A/i-pair and through self-energy contributions 
which contains the A. Our scheme puts special emphasis on the evaluation of Pauli exchange 
terms. 

By the employment of an effective Fermi momentum we have evaluated the response for 
40 C a at several momentum transfers, obtaining a good agreement with data. These results 
were obtained as a consequence of the addition of several ingredients, which basically can be 
grouped into two sets: the opening of new decay channels (as 2p2h, lplA2/i, etc.; represented 
by the second order self-energies) which reduce the height of the quasi-elastic peak and spread 
the response, but moves the peak towards lower energies. And in the second group we put 
the Hartre-Fock self-energy and the RPA, which moves the peak to higher energies, but do 
not provide the required reduction of the intensity. In both cases the action of the exchange 
terms and the election of the interaction is very important. 

In many works, it is mentioned that the longitudinal response seems more elusive than the 
transverse one, even though the underlying physics is more complex in the last case. In view 
of our results, we think that the just mentioned assertion should be revised. Referring this 
affirmation we recall that, i) the first experimental results which separates the longitudinal 
and transverse channels [2], underestimate the longitudinal response, ii) the free Fermi 
response strongly differs with data only in the longitudinal case. Hi) A very simple model, 
which uses the ring approximation with a ir + p plus the Landau-Migdal ^'-interaction gives 
a good adjustment for the position and height of the transverse peak. The conjunction of 
these elements certainly shows the transverse channel as more understandable. However, a 
word of caution should be given: this simple model provides no intensity in the dip region and 
when all possible mechanisms (mesons exchange currents, initial and final state interactions, 
etc.) are incorporated, the resulting response shows a tendency to overestimate data. Our 
model is more complex than the ring series with the g'+n + p-interaction. But it is simpler 
than the more elaborate ones for the transverse channel already mentioned. These points 
suggest that it is precisely the transverse channel the most difficult to be understood. 

As a final comment, we would like to stress that the many-body problem in inclusive 
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quasi-elastic electron scattering is quite involved. Any result is strongly dependent on the 
nuclear interaction, which is full of unknown coupling constants and parameters. Recent 
improvements (see for instance, Ref. [28]), show that we are closer to the solution of the 
puzzle, which necessarily implies a complex many-body calculation together with an ade- 
quate election of the nuclear interaction. In any case, the simultaneous study of several 
momentum transfers should be given. In this sense, we think that the present contribution 
is a step further in this direction. Clearly, our next objective is the inclusion of ground state 
correlations and meson exchange currents in order to analyze the transverse channel. 
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APPENDIX A 



In this appendix we give the explicit expression for C(q, cu, m p , rrih) and Ca(q, toa, rrih) It 
is convenient to change variables using the dimensionless quantities Q — q/kF, h — h/kF, 
v = u>/2e F , A p = m/m p , X h = m/m h , A A = m/m A , c = and 5 = fr^zj^; e F = 

k F /(2m) is the Fermi energy and m is the nucleon mass. 

C(Q, v, A p , X h ) = mk F J -^0(\h + Q\- 1)9(1 - \h\)5(u - (X p (h + q) 2 - X h h 2 )) (50) 



If \ p = Xh = 1, 



w^^h" (51) 



for Q < 2 and < v < Q - \Q 2 ; , 

„ „ x mkF 1 , ,v C?,o, 

C(Q,,,l,l) = ^-(l-(--^r) (52) 

for Q < 2 and Q - \Q 2 < v < Q + \Q 2 or Q > 2 and -Q + \Q 2 < v < Q + \Q 2 and 
£(Q, v, 1, 1) = otherwise. 
If A p 7^ X h , 

C(Q,is,X p ,X h ) = (53) 
when \X h - Ap + 2i/ - A p Q 2 |/(2 A P Q) > 1 or v < (X p - X h )/2, 

C(Q, is, A p , X h ) = ^- 2 ^(2u+ (X h - X p ))/(X h X p ) (54) 
for Q < 2 and < v < (X p - X h (Q - l) 2 )/2, 

C{Q,v,\,,\ h ) = — [{X h - Ap) 2 - Q 2 A p (A fe + Ap) + 2u(X h - X p ) + 

+2gA pV /A fe Apg 2 - 2z/(A fe - Ap)] \ (55) 

for Q > 2 or v > (X p - X h (Q - l) 2 )/2. 

Now, we give the expression for Ca(Q, v , Aa 5 A/J, 

r rl 3 h 

Ca(Q, v, A A , A,) = mkp J -^0(1 - \h\) (h 2 Q 2 - (h ■ Q) 2 ) 5(u - (c(h + q) 2 /2 - h 2 /2 + 5)) 

(56) 
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It is convenient to re-defined variables as, 



then, we get, 



if | (ai + 03)702! > 1 and 





1-c 


Oi = 


2 


a 2 = 


-c|Q 




c 


a 3 = 


z/ 



- * (57) 



£a(Q, v, Aa, A/j) = (58) 



2-7T 

£>a(Q, Aa, A/j) = — {-^/6 + (^-2a 1 a 3 )/4-c 2 /2 + 



+a^Ln/6 - (a| - 2o ia3 )Cn/4 + c 2 ^Ln/2} (59) 



if \(ai + a 3 )/a 2 | < 1, with h min = \(a 2 - v/ a l - Aa 1 a 3 )/(2a 1 )\. 
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APPENDIX B 



In this appendix we show the residual interaction V when at least one A is involved and we 
present the explicit expression for the Fock self-energy for the A, T, A (p). 

We have employed the following Viai/i,i/i'ia' interaction (this interaction is required in 
the evaluation of the fourth diagram in the first line of Fig. 2 and in the second one of the 
second line, already called as Q 2 ), 



VW,iA'i/,'(p) = %^ Kan(p) (9'an(p) <? S + h' AN (p) a ■ p S ■ p) r ■ T (60) 



where, 



9' an{p) — 9 an ~ 7ZS 



Kan(p) r P 2 



and 



hi i \ p2 T Ian(p) n P 2 fM v 

hMp) = ~fta + ^m c > an p^; (61) 

r / s ^TrAN,pAN ~ (^,p) 2 / fi0 N 

1 tt AN, p AN \P) = A 2 , 2 ' ^ bI > 

ir,p ' P 

Analogous expressions are used for ViAi/i,ip'i/i' (needed in the evaluation of the first diagram 
of the second line in Fig. 2, called Qi), ViAih,iA'ih' (Q3) an d ^iai/i,2A {Qa}- For each of these 
interactions, the coupling constants and parameters should be replace by their corresponding 
NN, NA and AA values. Also Pauli matrices must be replaced by the corresponding 
transitions matrices S and T or S and X, the 3/2-3/2 spin matrices (see ref. [29]); depending 
on the character of the mesonic vertex. 

For the parameters entering into our theory we have set at 140 MeV (770 MeV) the mass 
of the pion (rho meson). The pion coupling constant f^ NN /An=0.081, f^AN^ f-n-NN and 
fnAA= I fnNN- For the rho meson we have used C p nn = C p an = C p aa = 2.18. The different 
mesons cut-offs at the vertices are set to A wNN = 1300 MeV/c, A pNN = 1750 MeV/c while 
all the remaining cut-offs are set to 1000 MeV/c. For the Landau Migdal parameter g' AN 
and g' AA we have taken the values 0.55 and 0.40, respectively. 

We present now the expression for E^(p), 

s a(p) = o ^ {(3g f AN ~ l)vo(A^AJv,/i7r,p) - 2C pAN v (A pAN ,n p ,p) 

6 A*7T 

+v c (A nAN , p) + 2 C pAN v c (A pAN , ii p , p)} (63) 
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where we have introduced the functions v and v c , defined as, 

MA ,,, P) = /^ 9(1 _ W)( _^_ } . 

<K»p) = / <g* s(1 " 7 r +Wl>? < a4'(*-> (64) 

which can be easily integrated. 
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Table I. 



q [MeV/c] 




d(1) 






d(1) 
n Ah 






X p 


x h 


Alo [MeV] 


Aa 


x h 


Alo [MeV] 


300. 


1.33 


1.54 


11.0 


1.08 


1.11 


4.0 


325. 


1.32 


1.54 


11.0 


1.06 


1.10 


3.5 


350. 


1.31 


1.54 


9.0 


1.05 


1.09 


3.0 


400. 


1.28 


1.54 


9.0 


1.05 


1.08 


3.0 


450. 


1.32 


1.56 


9.0 


1.04 


1.07 


2.5 


500. 


1.37 


1.59 


7.0 


1.04 


1.07 


2.0 



Table I: Effective masses and energies shifts employed in R p ^(q,uj) and R Ah (q,Lo). We 
have used dimensionless quantities X„ H = m/rrii, where % = p, h or A. 
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Figure Captions 



Figure 1: Goldstone diagrams stemming from Eq. (42). In every diagram two wavy lines 
represent the external probe with momentum and energy q and cu, respectively. The 
dashed line represents the residual interaction. Continuous lines represent either a 
particle or a hole, while a double-continuous line represents the A. Only forward-going 
contributions are shown. 

Figure 2: The same as Fig. 1, but for Eq. (49). 

Figure 3: The R^h response at momentum transfer q = 400 MeV/c. The short-dashed 
line represents the free response of Eq. (47), the long-dashed line the R^ h response of 
Eq. (48), while the continuous line is the final result R^ h of Eq. (49). 

Figure 4: Different contributions to the R p h response. In each panel the short-dashed line 
represents the free response, the continuous line is the final result of Eq. (42) and the 
meaning of the long-dashed line for each panel is explained in the text. The momentum 
transfer is 400 MeV/c for all panels. Data was taken from Ref. [3]. 

Figure 5: The R p h response at different momentum transfer. Short-dashed line are the 
free responses and the continuous ones are the sum of R p h plus R^h- The momentum 
transfer in each panel is: (a) 300 MeV/c, (b) 325 MeV/c, (c) 350 MeV/c, (d) 400 
MeV/c, (e) 450 MeV/c and (f) 500 MeV/c. Data was taken from Ref. [3]. 
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Fig. 1 




Fig. 3 




Fig. 4 
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